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Abstract 

We study a discrete-time resource flow in Z d , where wealthier vertices attract 
the resources of their less rich neighbors. For any translation-invariant probability 
distribution of initial resource quantities, we prove that the flow at each vertex 
terminates after finitely many steps. This answers (a generalized version of) a 
question posed by van den Berg and Meester in 1991. The proof uses the mass- 
transport principle and extends to other graphs. 

AMS Subject Classifications: 60K35, 90B15, 68M14. 
Keywords: clustering, mass transport, cellular automaton. 



1 The model and results 

We consider the following model for distributed clustering. Initially each vertex x G Z d 
is assigned a random amount of resource, ^ C (x) ^ oo, sampled according to a 
translation-invariant distribution. We denote x ~ y if x, y are adjacent on Z d , write 
x ~ y if x = y or x ~ y, and take Q x = {y : y ~ x}. At step n — 0, 1, 2, . . . , each vertex 
x holding some amount of resource will transfer all its resource to the vertex a n (x) G Q x 
with the maximal amount of resource. All the vertices update simultaneously, leading 
to the state (C n+ i(x),x G where C n+ i(x) is the sum of resources transferred to x 
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at step n. More precisely, we take a n (x) = argmax{C ri (?/), y G Q x } if C n (x) > 0, and 
otherwise If there is more than one y G Q x that attains the maximum, we 

say there is a tie at x and in this case a n (x) is chosen uniformly at random among the 
maximizing vertices. Finally we take E n (y) = {x G Q y : a n (x) = y} and C n+ i(y) = 



Note that, except for the possible tie breaking during the dynamics, all the randomness 
is contained in the initial data. Let P and E denote the underlying probability measure 
and its expectation, for both the initial resource quantities and possible tie breaks. 

This model is a simple example of a self-organized structure that emerges from a dis- 
ordered initial state. Instances of this type of phenomenon in several fields of science 
are mentioned in [I], where the model was introduced. 

We are interested in the following phenomena, concerning the stability properties of 
this dynamics. 

Question 1 . Does each vertex transfer its resource eventually to the same fixed vertex? 

Question 2. Does the the flow at each vertex terminate after finitely many steps? 

Question 3. If the answer to the previous question is affirmative, is the expected value 
of the final resource equal to the expected value of the initial resource? 

Van den Berg and Meester [3] considered this model on Z 2 with continuously-distributed 
i.i.d. initial resource quantities and answered Question [TJ Namely, they showed that 
F[a n (x) is eventually constant] = 1. For the case of i.i.d. initial distributions supported 
on N (see [3] for the precise hypotheses), they also answered Question [2] in any dimen- 
sion, i.e., they proved that F[a n (x) = x eventually] = 1. 

Again for i.i.d. continuous initial distributions on the two-dimensional lattice, van den 
Berg and Ermakov pQ considered some percolative properties of the configuration after 
one step, and reduced Question [3] to a finite computation. As the calculation would be 
too big, even for the most powerful computers, they performed Monte Carlo simulations, 
obtaining overwhelming evidence that the answer to this question is positive. 

In this paper we answer Question [2] in a rather general setting. We consider any 
dimension and allow any translation-invariant initial distribution. Question [3] remains 
open. 

Theorem 1. On Z d , for any translation-invariant distribution of the initial resources, 
the flow at each vertex almost surely stops after finitely many steps. 

The proof of local fixation essentially consists of ruling out, one by one, all the possibil- 
ities which could lead to a different scenario. For this goal we use the mass-transport 
principle in different ways. 
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In Section [2] we prove Theorem [TJ In Section [3] we discuss generalizations of our results 
and conclude with a remark about Question [31 



2 Proof 



Theorem [T] will be proved using a combination of lemmas. We begin by introducing 
some extra notation. 

Let | • | denote the cardinality of a set. For each x G Z d and n = 0,1,2,..., one and 
only one of the following events will happen (see Figured]): A n (x) = \C n (x) = 0], 



B n (x) = [C n (x) >0,E n (x) = {x} 
E n (x)\ >l],£ n {x)= [E n (x) = & 

C' n (x) 
L n+ i(v) 
S n (w) 



, C n (x) = \E n (x) = {z} for some z ^ x\, V n {x) 
. For n — 0, 1, 2, . . . we also set 

C n (a n (x)), 

a n (L n (v)), with L (v) = v, 
{v : L n (v) = w}. 



Observe that C' n (x) ^ C n (x). L n {v) is the location at time n of the resource that 
initially started at vertex v, and S n (w) denotes the set of all vertices whose initial 
resource is located at vertex w at time n. 
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Figure 1: Schematic diagram for the notation used in this paper, omitting the index n. 



Lemma 1 ([3], p.337). For each x, F(V n (x) i.o.) = P(£„(x) i.o.) = 0. 



Proof. We repeat the proof of [3] for the convenience of the reader. 

First, E n (x) can happen for at most one value of n, after which one always has A n {x). By 
translation- invariance one has E( | E n (x) |) = J2 W F[a n (x+w) — x\ — ^2 w F[a n (x) = % — 
w\ = 1. This is the most common use of the mass transport principle. Defining D n (x) = 
\E n (x)\ - 1 gives E[D n (x)] = 0. So F[D n (x) > 0] < E[D n (x)+] = E[D n (x)~] = 
F[D n (x) = — l] . But the last event corresponds to £ n (x), therefore J2 n W[D n (%) > 
0] ^ 1. Since \D n (x) > 0] corresponds to V n (x), the result follows by the Borel- 
Cantelli Lemma. □ 

Note that C n+1 (x) > C n (x) can only happen on the event V n {x). So, almost 
surely, (C n (x)) is a nonnegative, eventually non-increasing sequence, and therefore 
lim^oo C n (x) exists. 

The following lemma permits us to extend the results of [3] to a general translation- 
invariant probability distribution for the initial resource quantities. 

Lemma 2. F-a.s., ties cannot happen infinitely often at a fixed vertex. 

Proof. For a pair of vertices y ~ x, let 

T n (x, y) — [0 < C n (x) = C n (z) = max {C n (tf), w ~ ?/} for some z ~ 2 7^ x\ 

denote the event that there is a tie at y at time n, and x is one of the candidates for 
its resource. Denote by F y>n the cr-field generated by 

(C m (z),m ^ n and z E Z d ) and (a m (,z),m < n and z E Z d or m = n and z^i/), 

i.e., J-y jn contains all the information up to step n, except for the possible tie breaking 
at vertex y at the n-th step. Write M^\ y = \E n (x)\{y}\. Note that MV and T n (x,y) 
are J^-measurable. Conditioning on J-^n, if MjL — and T n (x,y) occurs, then 
with probability at least ~, a n (y) 7^ x, in which case S n (x) will happen and thus 
C m {x) = for all m > n (therefore T n (x,y) can never happen again). If M?X > and 
T n {x, y) occurs, then with probability at least ^tj, a n (y) = x, in which case T> n (x) will 
happen. Therefore, the occurrence of T n (x,y) for infinitely many ra's implies almost 
surely the occurrence of T> n (x) infinitely often. But the latter event has probability 
by Lemma [TJ □ 

Lemma 3 ([3J, p. 338). F-a.s. for each vertex x, only one of the events limsup n ^4 n (x), 
limsup n B n (x) or lim sup n C n (x) will happen. 
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Proof. Once we know that there are finitely many ties at most, the argument of [3] can 
be applied. We present the proof for the convenience of the reader. 

By Lemmas [1] and [2] we can take no so that neither T> n (x) nor S n (x) will happen, nor 
will there be a tie at any y G Q x , for any n ^ n . If A n (x) happens for some n, it 
also happens for all m > n. It is thus enough to prove that, if C n (x) happens for 
some n ^ n , then C n+ i(x) also happens. Now, if C n (x) happens, we have C n+ i(z) ^ 
C n (x) > C n {y) = C n+ i(x), where E n (x) = {y} and z = a n (x). The strict inequality 
holds because there cannot be a tie at y. Therefore a n+ i(x) ^ x and, since V and S 
have been ruled out, we must have C n+ i(x) again. □ 

By Lemma [3] we can say that each vertex is uniquely either an A- vertex or a B- vertex 
or a C- vertex. 

Corollary 4. F-a.s., C' n (x) - C n (x) -> for all x G Z d . 

Proof. If x is an A- vertex or a B- vertex then C' n (x) — C n (x) = eventually. So suppose 
x is a C- vertex. Take Uq so that C n (x) happens for all for n ^ n . By Lemma [1] we can 
further assume that T> n (z) will not happen for any z G Q x and n ^ no, in particular 
C n+ i(z) ^ C n (z). We claim that C' n+2d (x) ^ C n (x) for all n ^ no- Since in addition 
C n+ 2d(x) ^ C' n+2d (x) and C n (x) converges, this will finish the proof. 

Let n ^ no be fixed and for m ^ n let V„ = {y ~ x : C m (y) > C n (x)j. Since 
resource quantities no longer increase, C V m . Moreover, the occurrence of C n (x) 

implies that a n (x) ^ x and, since P n (a n (x)) cannot occur, a n (x) is a C- vertex and 
E n (a n (x)) = {x}. Thus, C n+ i(a n (x)) = C n (x) and therefore a n (x) G" Ki+i- If Ki+i = 
we are done, and if V n+ i ^ we find again that a n+ i(x) G Ki+i\Ki+2- Proceeding this 
way, and since < 2d, we must find V n+ j = in less than 2d steps. □ 

Lemma 5. For each fixed k G N and x G Z d , 

P (0 < C' n {x) - C n (x) < 5, \S n (x)\ ^ k) -»■ 
as 5 — >■ 0, uniformly in n. 

Proof. For f G Z d consider S n (L n (v)), the set of vertices whose initial resource were 
joined with that of t> by time n (according to the dynamics rules, once two or more 
initially distinct sources join, they never separate), and define N n (v) := \S n (L n (v)) |. 
For each fixed v, N n (v) is non- decreasing in n. 

Let A n {y) = C n (L n (v)) be the total amount of resource at L n {y) at time n and let 
A' n (v) = C' n (L n (vy\ = C n (L n+ i(v)) . For each vertex v, the value A' n (v) — A n (v) is 
nonnegative and it can only decrease when N n (v) increases. 
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It is therefore the case that inf [^(f ) — , the infimum taken over all n such that 

N n (v) ^ k and A' n (v) > A n (v), is a strictly positive random variable. In particular, for 
fixed k, 

P(0 < A' n (v) - A n {v) < 5, N n (v) ^ k) ->■ as 5 -> uniformly in n. (1) 

We shall relate the above limit to the desired result via a quantitative use of the mass- 
transport principle, what we call an unlikelihood transfer argument. 

Let us omit 5, k from now on and denote the event in (JT]) by U n [v). For x G Z d , 
let V n (x) denote the event that U n (v) occurs for some v G S n (x). If U n {y) happens 
for any such v then it happens for all of them (because the values of A n , A' n and N n 
are constant within S n (x) and equal C n (x), C' n (x), and (^(x)!, respectively). Writing 
m(v,x) = tL n ( v )=x,u n (v) we have that tu n (v) = J2 x m ( v ^ x )- ^ follows from the trans- 
lation invariance of the process that E^^m(to,w + y) = JsL^2 y m(w — y,w), which 
means 

E XL Tn(w, x) = EJ2 V m(v, w) V w, (2) 

giving 

F[U n (v)]=E[\S n (x)\t Vn{x) ]^¥(V n (x)) Vv,x. (3) 

Now V n (x) is exactly the event considered in the statement of the lemma and since 
P[W n (V)] as 5 uniformly in n the result follows. □ 

Proof of Theorem^- By Lemma [3] it is enough to show that P(a n (x) ^ x) — > as 
n — > oo. For any choice of 5 > 0, if a n (x) ^ x then either C' n (x) = C n (x), which implies 
that there is a tie at x, or C' n (x) — C n (x) ^ 5 or < C' n (x) — C n (x) < 5. By Lemma [2] 
and Corollary H] the probabilities of the first two events tend to zero as n — > oo. 

For any choice of k G N, the last event can be split into two cases, according to whether 
1^(^)1 ^ k or not. By Lemma [5j the probability that this event happens in the first 
case tends to zero as 5 — > 0, uniformly in n. Now, with m(v,x) = t Ln ^ =x , ([2]) gives 
E(\S n {x)\) = lso¥(\S n {x)\>k) <l 

So, if we consider limsup n _ >00 P(a n (x) ^ x) , let 5 — > and then let k — >■ oo we get the 
desired limit. □ 

3 Concluding remarks 

We conclude this paper by discussing how Theorem [T] extends to more general settings 
and why a positive answer to Question [3] does not follow from the previous arguments. 
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Generalizations Our proof applies in other settings with much generality, as long as 
the mass transport principle ([2]) is true. It thus covers cases when the distribution of 
initial resources is invariant with respect to a transitive unimodular group of automor- 
phisms. Examples include Cayley graphs, regular trees, etc. The only change in the 
proof is to replace 2d by the graph degree. 

It also covers graphs that are locally finite and can be periodically embedded in M. d . 
Namely, one can consider graphs whose vertex set may be written as [J] x Z d , where 
[J] := {1, . . . , J}, and whose edge set is invariant under the mappings (j, x) i— >■ (j, x + y) 
for all i/eZf Here translation invariance is understood as the distribution of the initial 
resources being invariant under the above mappings. 

The changes in the proof for the above case are the following. In Lemma [TJ notice 
that E n (j,x) will happen at most once for each j G [J] and fixed x; so D n (x) := 
(Y?i\E n {j^)\) ~ J satisfies E[ D n {x) + ] = E[E n A^)1 < J, thus D n {x) > 
can occur for at most finitely many n's; finally T> n (j, x) corresponds to \E n (j, x) | — 1 > 0, 
which in turn implies that D n (x) > or that S n (j',x) occurs for some other j' G 
[J]. In Lemma [2] and Corollary H] replace 2d by max.,- deg(j, x). In Lemma [5] we take 
m(v,x) = J2j j> ^L n (j,v)=(j',x),u„(j,v) and the mass-transport principle ([2]) holds with the 
same proof. In the proof of Theorem [TJ we can write E( ^ . \S n (j, x)\) = J, giving 
F(\S n (j,x)\ > k) < ± for any (j,x). 

Open question It would be natural to try to answer Question [3] using an argument 
similar to the proof of Theorem [TJ in a way that uses the result to reinforce itself. 

First notice that a yes to Question[3]is equivalent to showing that every amount of initial 
resource eventually stops moving. Indeed, if C^x) := \im n C n (x) and F v denotes the 
event that L n {y) is eventually constant, then E[Coo(x)] = E[Co(f )1f v ] , so a yes to 
the question is equivalent to F v happening almost surely. This equality follows from 
Theorem [TJ and with 

/ \ I Cq(v), if L n (v) = x eventually, 
m(v , x) = < 

10, otherwise. 

Let us denote by U n (v) the event that, for some m > n, L m+ i(v) ^ L m {v). This event 
means that the resource that started at v is still going to move after time n. Therefore 
Question [3] boils down to whether we can show that P[Z4(i>)] -»■ as n — > oo. 

Let V n (x) denote the event that the resource found at vertex x at time n will leave x 
after time n, that is, a m (x) ^ x for some m > n. Note that this is equivalent to the 
event that U n {v) occurs for some v G S n (x), and again if it happens for any such v then 
it happens for all of them. This implies that the analogue of ([3]) holds. 
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By Theorem [TJ P[V n (a;)] — > 0. The inequality in the analogue of (^j does not help 
here, but we could make use of the equality if the distribution of 1^(^)1 satisfied an 
appropriate moment (or uniform integrability) condition. This is another quantitative 
use of the mass-transport principle, that we call unlikelihood sharing principle. 

For instance, suppose one can prove that E(\S n (x)\ a ) stays bounded as n — >■ oo for 
some a > 1. Then P[V n (x)l —> implies that P[ZY n (u)l — » since 

E(Slv) = P(V)E(S|V) ^ P(V)[E(,S a |V)] 1/a < [P(V)] 1_ -E[5 a ] 1/Q , 
where S = \S n (x)\ and V = V n (x). 

Unfortunately we do not know how to control the tail of the distribution of jS^a;)!, 
and Question [3] remains an open problem. 

While the final version of this paper was in preparation, it was found by van den Berg, 
Hilario and Holroyd [2] that Question [3] as stated has a negative answer: there do exist 
nonnegative translation-invariant initial distributions for which resources do escape to 
infinity. 
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